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Preface 

This  report  is  the  seventh  concerned  with  research  accomplished 
in  connection  with  Navy  Contract  Horn*- 433 (00),  between  Dunham  Laboratory, 
Yale  University,  and  the  Office  of  Naval  Research.,  Department  of  the 
Navy.  In  this  report  is  given  a discussion  of  the  solutions  for  a 
pair  of  simultaneous  nonlinear  differential  equations  that  may  apply 
to  phenomena  of  interest.  These  equations  are  studied  analytically, 
and  particular  examples  are  solved  with  an  analog  computer. 

The.  research  was  carried  or.  and  the  report  written  by  the 
undersigned. 

K.  J.  Cunningham 
New  Haven,  August  1954 
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Abstract 

A pair  of  simultaneous  nonlinear  equations 

dx/dt  = (<^AX)  [kx  - x ' fx(y)jx 
dy/dt  = ((Xy/k,)  |k  - y - f_(x)Jy 

may  represent  problems  of  in+erest  involving  certain  biological  or 
physical  phenomena.  These  equations,  together  with  several  special 
cases,  are  investigated  analytically  and  information  about  their 
solutions  is  obtained,  A variety  of  dif.'.'orent  solutions  can  occur, 
dependent  upon  the  coefficients  in  the  equations  and  upon  the  coupling 
functions  f and  f ...  Criteria  are  developed  from  which  the  properties 
of  the  solutions  can  be  predicted.  Several  numerical  examples  are 
solved,  making  use  of  an  electronic  analog  computer. 
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i. 

I.  Simultaneous  Growth  Equations 

A problem  which  has  attracted  the  attention  of  both  mathematicians 
and  biologists  is  that  of  describing  mathem#  ' ically  the  effects  of 
environment  upon  the  population  of  a species  of  animal.  One  phase  of 
this  problem  involves  the  effect  of  competition  between  two  different 
species.  Where  tne  two  species  ©Insist  on  a common  food  supply,  or 
where  one  species  is  eaten  by  the  other,  or  where  additional  effects 
exist  so  that  the  population  of  one  species  influences  that  of  the 
other,  a mathematical  description  of  the  situation  must  require 
simultaneous  differential  equations . In  general,  these  equations  are 
nonlinear  and  are  more  or  less  complicated  depending  upon  the  number 
of  effects  that  are  considered - 

Volterra'1'  has  studied  in  scane  det^-il  the  case  where  the  competition 
equations  are 

dx/dt  = x » (aykv)(kx  - p^x 

dy/dt  “ y = (a/kJOq,  - fi  x)y 
j j ./  j 

Here-  t is  the  independent  variable  time,  x and  y are  dependent 

variables  representing  the  two  populations,  and  k , a , p , k , a , B 

x x x y y y 

are  real  constants.  Coupling  between  the  variables  comes  about  through 
the  nonlinear  product  term,  xy,  in  each  equation.  The  variables  can 
bs  separated  in  these  equations  and  a solution  found  by  a process  par liy 
analytical  and  partly  graphical.  The  nature  of  the  solution  depends 
upon  the  coefficients  in  the  equations.  Among  the  possibilities  are  that 
one  species  disappears  leaving  th-_  other,  or  that  oscillations  occur 
in  the  populations. 

1.  V,  Volterra,  La  Lutte  pour  la  Vie,  (Gauthier-Villars,  Paris,  1931), 
ch.  X 
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Somewhat  more  complicated  equations  have  been  investigated  qy 


X = (a3AxKkx  - X “ 

y = («/V(ky  " 7 “ V)y* 

The  extra  terms  in  these  equations  provide  a sort  of  damping  effect  and 
insure  that  neither  population  goal  to  infinity.  The  method  of  analysis 
used  by  Volterra  is  no  longer  applicable.  Again,  a variety  of  solutions 
is  possible,  depending  upon  the  coefficients. 

A still  more  complicated,  pair  of  equations  has  been  suggested  by 


Hutchins onJ  as 


i = W&x'**  fX(y)JX 

v * (a  /k  )lk  - v - ? (x)lv 

* v » tr  • v _j- 


where  f (y)  and  f (x)  are  the  functions  coupling  one  variable  to 
x y 

the  other.  These  functions  are  generally  well-behaved  mathematically , 

and  it  is  required  that  f (0)  ^ 0 and  f (0)  = 0.  Tho  functions 

x y 

might  be  polynomials  of  the  sort 

fx(y)  " pxy  + + s^y3 

f (:c)  « p x + y x2  + 8 x^. 
y y y y 

The  equations  of  Gause  and  Uitt  contain  only  the  first  terms  of  these 
polynomials  * 

Equations  of  this  same  general  sort  might  arise  in  describing 
other  type s of  physical  phenomena,  For  example,  certain  1 rinds  of 
chemical  reactions  progress  at  a rata  that  depends  upon  the  amount  of 
each  component  present  in  the  reaction.  Simultaneous  equations  describing 
the  amount  cf  each  component  are  quite  similar  to  the  competition 
equations.  Similarly,  an  electrical  system  can  be  conceived  which  also 


2,  G,  F.  Gause  and  A.  A.  Witt,  American  Naturalist,  69,  596,  (1935) 

3,  G.  E.  Hutchinson,  Ecology,  2S.  31b.  (1947) 
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would  be  described  by  the  equations.  The  two  voltages  of  a pair  of  d-c 
generators  sight  build  up  in  this  way.  if  the  field  windings  were 
suitably  connected.  Two  field  windings  on  each  inachine  would  be 
required,  one  excited  by  the  machine  itself  and  the  other  excited  from 
a cross  connection  with  the  seconu  machine* 

In  general,  the  equations  apply  to  any  situation  where  two  effects 
tend  to  grow  with  time,  but  the  rate  of  growth  of  each  is  influenced 
in  some  way  by  the  other. 


II,  Most  General  Form  of  Equations 
11,1  Analysis  of  equations 

The  pair  of  simultaneous  equations  that  arc  considered  here  of 
the  most  general  form 

dx/dt  - x = (a^/k..)  - x - f_Xyj]x  (1) 

dy/dt  - y « (cLy/ky)  jlcy  - y - fy(x)Jy,  (2) 

In  these  cq-dtions,  t is  the  independent  variable  and  usually  represents 
time,  x and  y are  the  two  dependent  variables,  and  &x,  ay,  kx,  ky 
are  real  constants.  Functions  fx(y)  and  fy(x)  are  continuous, 
single-valued,  well-behaved  functions  that  can  be  differentiated  with 
respect  to  their  arguments.  Both  functions  vanish  for  zero  argument, 
fy(0)  = 0 and  f (0)  = 0,  In  many  cases  the  functions  take  the  form 
of  polynomials  such  as 

Uy)  ° Pxr  + (3) 

fy(x)  = 3yx  + yyX2  (4) 

where  B , 0 , y , y are  real  constants.  It  is  through  those 
x y x y 

functions  that  coupling  exists  between  the  two  dependent  variables# 


\ 
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If  there  is  no  coupling  between  variables  x and  y,  f (y)  a 0, 
and  Eq.  (1)  becomes  the  Verhulst-Pearl  equation^ 

x = - x)x.  (5) 

This  is  an  example  of  a Bernoulli  equation  and  has  the  exact  solution 

x - [k^"1  + (x^1  - k^”1)  expC-q^t)]'1  (6) 

where  x ■ x at  t *»  0,  Curves  for  x as  a function  of  x and  for 
o 

x as  a function  of  t are  shown  in  Fig.  1*  with  several  initial 
conditions.  If  |xp|«|lqcj,  approximately  x = xq  exp(a^t)  and  the 
solution  starts  off  as  it  v/ould  for  a simple  growth  equation  having  a 
constant  growth  factor  a _.  With  a positive>  x always  approaches 

X X 


iroiiio  ^ rsQwrrd}  ftce  v/h0  whar  v 


r»  r-\r\  o Td  r»r»  nnf'i 


If  xq  is  negative,  x pauses  through  infinity  and  becomes  positive 
so  as  to  approach  a positive  k^.  The  value  x = 0 is  thus  a point  of 
unstable  equilibrium,  while  the  value  x = kx  is  a point  of  stable 
equilibrium. 

With  no  coupling,  f (x)  a 0,  also,  and  exactly  the  same  sort  of 

y 

solution  applies  for  the  equation  in  y as  has  just  been  discussed  for 
the  equation  in  x. 

In  the  more  general  case  where  coupling  between  the  variables  is 
present,  it  is  necessary  to  consider  the  complete  form  of  Eqs.  (l)  and 
(2).  The  nature  of  the  possible  solutions  for  these  equations  is  most 
easily  studied  by  considering  the  single  equation  obtained  as  their 
ratio, 

dir  (vV  - y - fyU)]y 

* ' WR,  - * - ' (7) 

The  independent  variable  t has  disappeared  in  writing  this  equation. 


4,  A,  J.  Lotka,  Elements  of  Physical  Biology,  (Williams  and  Wilkins, 
Baltimore,  1925).  p,  64 


r 

There  are  certain  singular  points5  for  Eq.  (7)  corresponding  to 

points  of  equilibrium,  where  both  numerator  and  denominator  vanish 

simultaneously.  The  values  of  x and  y at  one  of  these  singularities 

are  designated  as  X ana  y * The  singularities  can  be;  separated  intn 

s s 

four  groups  as  follows. 

Gr.  1,  xs  = 0,  ys  - 0 

Gr*  2*  xa  " V ys  - 0 

Gr,  3.  xs  =0,  ys  = ky 

Gr-  4-  x„  ■ ‘Sc  “ fx(y)-  y8  ■ 'v  " fy(x) 

The  first  three  of  these  singularities  are  reminiscent  of  those  which 
would  occur  for  a pair  of  Verhulst- Pearl  equations  with  no  coupling. 

The  fourth  group  may  contain  none,  or  a number  of  singularities, 
depending  upon  the  properties  of  functions  f and  f . 

x y 

The  nature  of  solutions  for  Eq.  (7)  near  each  singular  point  can 

be  explored  by  replacing  x with  (x  + u)  and  y with  (y  + v), 

s s 

where  u and  v are  snail  changes.  This  substitution  gives  the 
equation 

^ = V^Vs  + kyV  ~ ys‘]  - 2ysV  ~ ysfy(zs)  " ^'sV  (xs)u  ~ fyUs)v] 

- t ■"*  ' ft  . i 2 n..  ..  „ _ [ \ _ ■ t \ ' \ n 

QU  1(1  /K  / | K.  A T R u - a - u - a j.  v.y  ; • A i.'-  \V  /V  X V.V 

X X1-  X S X S S S X S S X s X S .3 


In  writing  this  equation,  series  expansions  for  the  coupling  functions 


have  been  used. 


fx(ys  + v)  « fx(yE)  + fx'(y8)v  + 
fy(xg  + u)  = fy(x8)  + fy»(xs)u  + 


v(ys}  = ^txM 

V (xs)  3 d/dx[fy(xs)]  • 


5.  N.  Minorsky,  Nonlinear  Hechanics,  (J.  W,  Edwards,  Ann  Arbor,  1947), 
Part  I, 
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Only  linear  terms  in  u arid  v have  been  retained.  At  any  singular 
point,  certain  terms  vanish. 

Hi  - xa  - fxM1.  ■ 0 


fie  - y - f (x  j*ly  = 0 
Ly  y s JJs 

so  that  Eq.  (3)  becomes 

*■  <vV  { bvV - * Dv ' } 

du  (ax/kx)  { [kx  - ^s  - Wl  U + L~xsfx'  (ys}]V  } 

. (?) 

This  equation  is  of  the  form 

dv  M Au  + Dv 
du  B Cu  + iv 

do) 

where 

A - -(«/ky)ysfy.(xs) 

(11) 

(12) 

0 ■ bAc’&x  - =*0  - fx(y3] 

(13) 

D“-<WVx'<^>- 

(U) 

The  nature  of  solutions  for  Eq,  (9),  .. nd  also  of  solutions 

near  the 

singularities  of  Eq,  (7),  depends  upon  the  characteristic 

exponents 

(\-  » 0/2)  [(B  * C)  ± UV2"| 

(15) 

whex-e  U “ (B  ▼ C;  ' •*-  4(aD  - BC),  it  is  necessary  to  examine  these 
exponents  near  e?oh  of  the  singularities. 


Gr,  1,  xs  - 0,  yg  » 0 
For  this  singularity 

A a o C « a 

x 

B - i D = 0 

y 

(B  + C)  *»  (a  + a ) 
x y' 

(AD  - BC)  **  -a  a 
x y 

U = (B  - C)2  > 0 


I 


8. 


i 


l 

i 

! 


If  a a < 0,  a saddle  point  exists;  if  a a > 0.  a nodal  point 
x y x y 

exists,  stable  if  both  ax  and  are  negative  and  unstable  if  both 

c and  a are  positive*  The  situation  is  summarized  in  the  stability 

•w  y 

diagram  of  Fig.  2. 

Gr.  2.  xs  - kx>  ys  =^0 

For  this  singularity 
A = 0 

B “ (cy  ky)Y  where  Y = ky  - fy(k,;) 

C a^x 
D " (0) 

(B  + C)  = (oY/l^  - ox) 

(AD  - BC)  - cxayY,rky 
U *»  (B  - C)2  > 0 

In  this  case,  the  quantity  Y is  important.  It  can  be  found 
easily  from  a graphical  plo*  relating  to  Eq.  (?)•  In  Fig.  3 is  plotted 
tha  curve  y = k - f .(*)  upon  axes  of  x and  y.  The  curve  is  the 

y y 

locus  of  those  values  of  x and  y which  mice  the  bracket  in  the 
numerator  of  Eq.  (7)  vanish.  Any  curve  representing  a solution  for 
Eq.  (7)  must  cross  this  locus  with  horizontal  slope.  Thus,  this  curve 
is  the  isocline  for  zero  slope,,  dy/dx  = 0,  or  for  y ° 0 in  Eq.  (2). 
The  ordinate  for  this  curve,  evaluated  at  x = k,.,  is  the  quantity  Y, 
and  iB  positive  in  Fig,  3» 

If  aya^Y/ky  > 0,  a saddle  point  exists;  if  axayY/ky<  0,  a 
nodal  point  exists,  stable  if  (B  + C)  < 0 and  unstable  if  (B  + C)  ”>  0. 
The  situation  is  summarized  in  the  stability  diagram  of  Fig.  h» 


L 


Fig.  2 Stability  of  Gr,  1,  Eqs  (1-2) 


fig.  5 Determination  of 


Fig.  6 Stability  of  'los.  (1-2),  Gr.3 
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Gr.  3-  x =0,  y **  k 

For  this  singularity 

A » -a  f » (0) 

y y 

B = 

y 

C - (a^/yx  where  X » kx  - fx(k.) 
D = 0 


(B  + C)  - (a//kx  - ay) 

(AD  - BC)  = axayX\ 

U - (B  - C)2  ^ 0 

This  case  is  analogous  to  the  one  just  preceding.  Here,  a plot  of 

the  curve  x = kx  - fx(y)  as  shown  in  Fig.  5 is  useful.  This  curve  is 

the  isocline  for  infinite  slope  of  a solution  curve  for  Eq,  (7), 
dy/dx  ■ oo,  or  for  x = 0 in  Eq.  (l).  The  abscissa  for  this  curve, 

evaluated  at  y * k^.,  is  the  quantity  X,  and  is  positive  in  Fig.  5. 

The  situation  is  summarized  in  the  stability  diagram  of  Fig,  6, 


3r-  *S  - ‘Sc  - ys_  “ rv(x5> 

The  number  of  singularities  determined  by  thooc  relations  depends 

upon  the  functions  f and  f , If  these  functions  are  nonlinear,  as 

for  example  Eqs,  (3)  and  (4),  numerical  determination  of  x and  y 

s s 

may  be  a fairly  tedious  process,  requiring  the  solution  for  the  roots 
of  an  equation  of  high  degree.  Probably  a simpler  and  more  informative 
approach  is  to  resort  to  another  grs  -cal  construction. 

Singular  points  In  general  are  located  at  the  intersection  of 
isocline  curves  corresponding  to  different  slopes,  i.e,  to  different 
values  of  dy/dx.  Thus,  if  the  isoclines  for  zero  slope,  y = 0 or 
y 3 k - i?  (x),  and  the  isoclines  for  infinite  slope,  x = 0 or 

y y 


I 


! 


f 

I 

l 


ii* 

x = - fx(y) , are  plotted,  singularities  are  found  at  their  inter- 

sections. The  number  and  location  of  the  intersections  vri.ll  depend 
upon  f and  f . 

A y 

Quadratic  functions,  such  as  Eqs,  (3)  and  (4),  lead  to  parabolic 
curves  for  the  isoclines,  as  in  the  example  of  Fig,  7,  This  figure  is 
like  a combination  of  Figs,  3 and  5,  but  here  the  functions  are  such 
that  both  X and  Y are  negative.  A maximum  of  four  intersections, 
and  thus  four  singularities,  can  occur  with  these  parabolas.  If  the 
coefficients  of  Eqs,  (3)  and  (4)  are  allowed  to  vary,  the  parabolic 
curves  change  in  shape  and  location,  and  the  number  of  singularities 
may  be  any  integer  from  zero  to  four  inclusive.  More  complicated  forms 
of  functions  f and  f might  lead  to  sven  more  singularities, 

J 

It  is  worth  noting  that  in  Fig,  7,  the  horizontal  axis  corresponds 
to  an  isocline  where  y «=  0,  and  the  vertical  axis  corresponds  to  an 
isocline  where  x * 0,  The  intersections  of  the  axes  vdth  each  other 
and  with  the  parabolic  curves  lead  to  the  first  three  singularities. 

For  the  singularities  of  the  fourth  group, 

a - -<y y y„v  <y  - -3I 
b - -cyyya  - -s 

0 - -<vV*s  ■ -R 
d - -<yy - -®1 

where 


R * vsA 
s s V/k: 


G = R/S  = a k x /a  kv 
x y s'  y xr  e 


m . <y(xs)  - Vdx[fy(xsj] 

m s y(y8)  - o/dyfqyj] 


Quantities  F.  and  S can  be  found  from  known  values  of  a ,,  a , 

x y 

k and  k together  vrith  values  of  x and  y taken  from  the 
x'  y s s 
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Fig.  8 Stability  of  Eos.  (1-2),  Gr 
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graphical  construction#  Quantities  M and  N can  be  found  from  the 

slopes  of  the  isocline  curves  for  y ■ 0 and  x.  = 0,  o.s  measured  at 

the  singularity.  The  isocline  for  y ■ 0 is  y =*  - f^.(x) , which 

has  the  slope  dy/dx  = — d/ dx[fy(xfj  » -f^1  (x) . This  slope  evaluated 

at  the  singularity  is  dy/dx 1,  = ~f  ’ (x  ) = -II  . Similarly, 

J(x  = x8;  y 8 

the  isocline  f or  x * 0 is  x = k - f r(y)  with  the  slope 

X X 

dy/dx  =*  ( do</ dy)"1  - |-d/dy  [fx(y)j  j'1  " f-f^Cy)]"1.  This  slope, 
evaluated  at  the  singularity  is  dy/dx],  = f-f  1 (y)l  ^ = -l/N. 

-Ky  * ys)  L '•  QJ 

Thu*.  M and  N can  be  determined  from  direct  measurement  of  the 

slopes  of  the  curves,  or  by  numerical  substitution  in  the  derivatives 

of  the  isocline  curves,  M = f 1 (x  ) and  N a f »(y  ). 

* y 's'  x w o' 

Important  combinations  of  the  coefficients  are  the  folio; ring. 

(B  + C)  « ~(R  + S)  =«  -S(G  1) 

(AD  - BC)  « Ii5(KN  - 1)  = S^G(MK  - 1) 

U - (H  - S)2  + 4RSMN  - 4S2g|]:g  + I/O) A - 1/2  + MrfJ 
The  following;  solutions  may  exist, 
a.  Saddle 


G(I4N  - 1)  > 0 


b.  Node 


Both  G(MN  - 1)  < 0 
and  g[(C-  + l/&)/4  - 1/2  ♦ MHj>0 
Stable  if  S(G  + 1)  > 0 
Unstable  if  S(G  + 1)  < 0 


C.  i'OCUS 


G[(G  + l/G)/4  - 1/2  + Mn)<  0 
Stable  if  S(G  + 1)  > 0 
Unstable  if  S(G  + 1)  < 0 

The  situation  is  summarized  in  the  stability  diagram  of  Fig.  G. 
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11,2  Numerical  examples  from  computer 
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in  order  to  check  these  conclusions  concerning  the  solutions  for 
Eqs.  (l)  and  (2)  a particular  pair  of  equations  war.  studied  with  an 
electronic  analog  computer.  The  equations  used  were  Eqs,  (l)  and  (2) 
with  the  quadratic  forms  of  Eqs,  (3)  and  (4),  These  quadratics  are  the 
most  complicated  functions  that  can  be  handled  with  the  available 
computer.  The  connections  for  the  computer  are  shown  in  Fig,  9«  The 
numerical  values  for  the  coefficients  of  the  equations  were 
a varied 

■Sc  - 3 
3X-  V2 

y = 1/6 

These  coefficients  yield  parabolic  isoclines  that  intersect  at  four 
points,  one  in  each  quadrant,  and  thus  give  a total  of  seven  singularities 
The  solution  curves  for  the  equations,  as  plotted  directly  with 
the  analog  computer,  are  shown  in  Fig,  10,  11,  and  12,  for  which 
ax  » +1,  -1,  and  -l/2,  respectively.  The  isoclines  for  y = 0 and 

i 

x - 0 are  shown  in  Fig,  10,  These  same  isoclines  apply  also  to 
Figs.  11  and  12,  since  a change  in  does  not  change  the  isocline 

curves.  Various  initial  values  of  x and  y were  used  in  each  case, 
and  the  resulting  solution  curves  are  shown.  The  direction  of  increasing 
time  is  indicated  by  the  arrowheads  on  these  curves. 

Important  numerical  data  applying  to  the  seven  singularities  for 
each  of  the  three  values  of  ax  are  given  in  Table  I,  These  data 
used  with  the  stability  diagrams  of  Figs,  2,  k,  6,  and  3 allow  the 
prediction  of  the  kinds  of  solutions  that  apply  near  each  singularity. 

The  predictions  are  listed  in  Table  I,  The  data  for  the  Gr.  4 singu- 
larities are  plotted  in  Fig,  13,  which  is  similar  to  Fig,  8* 
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Table  I 

Location  and  Types  of  Singularities  for  Examples 
x = (a2/kx)(kx  - x - - Y^y2)* 

y - (^y/ky)(ky  - y - &y*  - yyx2)y 

a ol 


a listed  below 
x 

k = 3 
x 


Bx-  V2 

'X  ' V6 


y 

py 


y - 

' V 


4 

1 

1/4 


1 1 

{ H 

i K 

srirgplAr 

point 

a 

V 

X 

s 

ys 

X 

Y 

M 

N 

i-2-y 

ft  i 

-j 

s 

G 

type 

j 

l.a 

1 

0 

0 

_ 

] 

!wsfc  node 

2. a 

] 

3 

0 

-1.25 

1 

| 

st  node 

i 

j 

3. a 

1 

0 

4 

-1.67 

j 

st  node 

i 

4:1  ;a 

1 

1.60 

1.76 

1.80 

1.09 

1.96 

*53  j 

.44 

1.21 

saddle 

4.2.  a 

1 

3.32 

-2.08 

2,16 

-.19 

-.41 

l.ll 

-.52 

-2.14 

saddle 

1 

4.3. a 

1 

-3.06 

4,72 

-.53 

2.07 

-1.10 

-1.021 

1.13 

-.86 

saddle 

! 

4 »4*a 

1 

-9.86 

-10.4 

-3.93 

-2.97 

11.7 

-3.29 

— • 61 

1.26 

saddle 

4 

! 

l.b 

-1 

0 

0 

saddle 

i 

I $ 

2.b 

-1 

3 

0 

-1.25 

saddle 

■ 

3.b 

-1 

0 

4 

-1.67 

saddle 

i n 

4.1. b 

-1 

1.60 

1.76 

1.80 

1.09 

1.96 

-.53 

•44 

-1.21 

unst  focus 

i « 

4.2.b 

-1 

3.32 

-2.08 

2.16 

-.19 

-.41 

-1.11 

-.52 

2.14 

unst  focus 

4.3. b 

-1 

-3.06 

4.72 

-.53 

2.07 

-1.10 

1.02 

1.18 

.36 

st  focus 

if 

4.4.b 

-1 

-9.86 

-10.4 

-3.93 

-2.97 

11.7 

3.29 

-2.61 

-1.26 

st  focus 

1 1 

l.c 

-v?. 

0 

0 

, saddle 

> W 

2.c 

-V2 

3 

0 

-1.25 

saddle 

’ 1 

3 . c 

-1/2 

0 

4 

-1.67 

saddle  j 

g* 

4.1. c 

-1/2 

1,60 

1.76 

1.80 

1.09 

1.96 

-.26 

• 44 

-.60 

st  focus 

i.  i 

4#2.c 

-1/2 

3.32 

-2.08 

2.16 

-.19 

-.41 

-.55 

-.52 

1.07 

uist  focus 

4.3 : c 

-1/2 

-3.06 

4.72 

-.53 

2.07 

-1.10 

.51 

1.18 

.43 

st  focus 

i | 

1 

• 

4 • 4 » c 

-1/2 

-9.36 

-10.4 

-3.93 

-2.97 

11.7 

1.6/, 

-2.61 

-.63 

unst  fccus 

st 


Fig.  13  Stability  of  Gr.  4 for  examples 
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The  types  of  solutions  predicted  in  this  way  are  the  same  as 
those  found  with  the  computer* 


II. 3 Sketch  of  solution  curves 

j. t is  possible  to  sketch  the  general  shapes  of  the  solution  curves 
without  having  recourse  to  a computer.  This  can  be  done  from  a 
knowledge  of  the  location  and  nature  of  the  singularities,  of  the 
isocline  curves  for  y = 0 and  x = 0,  and  of  the  asymptotes  for 
solution  curves  near  a node  or  a saddle.  Information  relating  to  the 
singularities  and  the  isoclines  has  been  given  above.  The  asymptotes 
are  considered  here. 

It  can  be  shovrn  that,  the  solution  curves  near  the  singularity 
corresponding  to  a node  or  a saddle  approach  a definite  slope  as  time 
approaches  either  plus  or  minus  infinity.^ As  tine  approaches  pin.? 
infinity,  the  slope  is 

dy/dxi  = m,  = A/(\^  - B)  = (A.  - C)/D.  (16) 

Jt  }+o? 

As  time  approaches  minus  infinity,  the  slope  is 

dy/dxj  " m2  = A/^2  " B)  ” ” G)D*  C1?) 

1-5*  — OD 

In  these  equations,  A,  3,  C,  D are  the  coefficients  from  Eq.  (10), 
is  the  more  positive  characteristic  exponent  and  is  the  less 

positive  characteristic  exponent,  as  found  from  Eq,  (15) . 

As  an  example,  a sketch  of  the  solution  curves  corresponding  to 
the  case  of  Fig.  10  is  shown  in  Fig.  14.  Isocline  curves  for  y = 0 

and  x = 0 are  plotted  free'!',  the  equations 

t 2 

y * 0:  y = k - p X - y X : y - 0 

«y  «)  v 

X = 0:  x = kx  - P^y  - y^y2;  x = 0 

6.  A.  A.  Andronow  and  C.  E.  Chaikin,  Theory  of  Oscillations,  (Princeton 
Uiiiv-riity  Pros?,  p^nn.Rt.on.  1949),  ch.  V. 

7.  B,  G„  Farley,  Proc.  IEE,  1+0,  1497,  (1952) 


! 
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Their  intersections  locate  the  seven  singularities,  the  natures  of 
which  have  been  oredicted.  in  Table  I. 

C: 

The  slope  of  a solution  curve  must  change  algebraic  sign  whenever 
it  crosses  one  of  the  isoclines,  y = 0 or  x ■=  0.  Thus,  with  the 
knowledge  that  the  origin  is  an  unstable  node,  so  that  the  solution 
curves  near  the  origin  in  the  first  quadrant  must  have  positive  slope, 
the  signs  of  the  slopes  of  the  solution  curves  in  each  region  can  be 
assigned,  as  shown.  Because  the  axes  are  also  particular  isoclines,  a 
solution  curve  can  never  cross  either  axis,  and  must  remain  within  the 
quadrant  in  which  it  starts. 

The  values  for  the  slopes  of  the  asymptotes  m^  and  are 

given  in  Table  II.  The  equality  of  a and  a in  this  example 

x y 

causes  m^  and  for  the  node  at  the  origin  to  be  indeterriinat0 . 

The  asymptotes  are  plotted  for  each  singularity  in  Fig,  14. 

The  asymptotes  are  aloo  separatrix  curves  near  a singulai’ity. 

Thus,  the  asymptote  nv,  near  Singularity  4.1  separates  those  solution 
curves  that  tend  toward  Singularity  2 from  those  that  tend  toward 
Singularity  3. 

k number  of  solution  curves  are  sketched  in  Fig.  14,  maidng  use  of 
the  information  collected  there.  These  solution  curves  represei.k  quite 
accurately  the  true  solutions,  as  found  with  the  computer  and  shown  in 
Fig.  10. 

A solution  curve  near  a focal  point  spirals  about  that  point  and 
approaches  it  from  no  definite  direction.  Thus,  there  are  no  asymptotes 
associated  -with  a focal  point,  and  this  aid  to  sketching  solution  curves 
is  not  available.  This  situation  applies  near  the  focal  points  of  the 

^ examples  shown  in  Figs,  11  and  12. 

1, 
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Table  II 

Asymptotic  Slopes  for  Example 
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III,  Special  Case  with  no  Squared  Terms 

The  form  used  in  writing  Eqs,  (1)  and  (2)  is  convenient  in  that  the 

parameters  k and  k represent  the  ultimate  values  of  x and  y 
x y 

when  there  is  no  coupling  between  the  variables.  This  form  is  not  well 


2 2 

adapted  to  apply  to  the  case  in  which  the  terms  in  x and  y"  are 


missing  on  the  right  sides  of  the  equations.  It  is  simpler  to  consider 
this  latter  case  by  itself,  rather  than  to  attempt  a modification  of 
the  preceding  analysis  to  include  it. 

The  equations  with  no  squared  terms  are 

* * (vVtjsc  - fx(y)]*  !l8) 

* “ WliV  ' Vx)]y  (19> 


where,  again,  f (0)  ■ 0,  f„(0)  = 0,  and  these  functions  can  be 

j 


differentiated.  The  ratio  of  these  two  equations  is 

1 . (VVCkv  ~ fy(xBy 

'‘VVGSc'  fx(y>> 


(20) 


As  before,  near  a singularity,  x is  replaced  by  (x  + u)  and  y by 

5 


(ys  f v),  giving 

* . <VV  { Ey5i:v'(?-6)]U  * LkV  * fV(Xs>1V} 

du  A>  { L'Sc  - W] u * !:  W ) 

where  fy;  (xg)  =*  d/dxjfy(xs)l  and  fx>  (yB)  **  d/<\v  rfx(ys)},  and  only 
linear  terms  are  retained.  The  coefficients  are 


(21) 


A = -(a  /k  )y  f 1 (x  ) 
y 7 s y s/ 


B » (a  J k ) 

y y 

C - (°/kx> 

D * -<“AKV  (v 

There  are  only  two  groups  of  singularities. 


fk  - f (x  n 

L y y s J 

A • fx(5'=0 


(22) 

(23) 

(24) 

(25) 


; 

3 

l 


26. 


Gr. 


xs  “ °'  ys 


«*> 


For  this  singularity 


The  stability  diagram  is  shown  in  Fig,  15- 


Or.  2.  y y„)  ■ kx,  fY(x8)  - kT 

For  this  singularity 

A - -(ay/ky)ygfy«  (xo)  C - 0 
B » 0 

Quantities  A and  D must  be  calculated  x'rom  the  equations.  The 

stability  diagram  is  shown  in  Fig-  16, 

An  example  of  this  special  case  is  that  for  which 

2 


5 - -(tAx)V,'(ya> 


fX(y)  = *2?  + V 

fy(x)  - Pyx  + yyx‘ 


(3) 

2 (4) 

2 


the  same  forms  used  previously.  If  the  quantities  y^K^/p^  and 
2 

y k/p  are  both  positi/e,  there  are  four  singularities  in  Gr,  2, 

J J'  v 

one  occurring  in  each  of  the  four  quadrants.  Furthermore,  the  algebraic 
signs  of  xg  and  fy*(xs)  are  the  same,  as  are  the  signs  of  ys  and 
f-.Ky,,)*  With  these  conditions,  the  sign  of  product  AD  is  the  same 
as  that  of  product  ayfiy • Typical  solution  curves  are  as  shovai  in 
Fig,  17,  With  other  relations  of  the  parameters  in  the  equations,  the 
solutions  ma.y  be  different,  of  course. 

If  in  Eqs,  (3)  and  (4),  both  - 0 and  yy  = 0,  the  original 
equations  reduce  to 


c 


4 " '“jAP&x  - vlx 

} - wCv  • vV 


(26) 

(27) 


t 


< 

\ 


Oc' 

■?,<  j • 

These  equations  are  those  studied  by  Volterra,  Only  two  singularities 
exist,  one  at  the  origin  and  the  other  at  the  point  xg  = k^/p^, 

ys  " VPx* 

IV.  Special  Case  with  one  Squared  Term 

A special  case,  intermediate  in  complexity  between  the  moot 


There  are  three  groups  of  singularities. 


Gr*  1*  x8  " °*  ys  0 


A = 0 


B = CL 


D = 0 


The  stability  diagram,  is  shown  in  Fig.  18, 


Gr.  2. 


k , y =0 
x*  •'  s 


-a 


D “ -Vx’ (0) 


A = 0 

B ‘ WfV  - 'Al! 

It  is  necessary  to  calculate  the  value  cf  B from  the  equation.  The 
stability  diagram  is  shown  in  Fig,  19. 


general  case  considered  first,  and  that  just  discussed,  is  the 

following. 

\ 

* - K/VOsc  - * - 

(28) 

2 

1 

1 

1 

* - WCV  - Vx3y  • 

(29) 

2 

The  x term  is  present  in  the  first  equation  but  there  is  no 

y2 

term  in  the  second  equation. 

l 

The  coefficients  applying  near  a singularity  arc 

• 

•••w.v11.1 

(30) 

\ 

b - vvOv  - yy] 

(31) 

4 

:i 

° ■ (vyfsc  - - yy] 

(32) 

i 

i 

i 

(33) 

I 
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29. 
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stable 
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saddle  *• 

Fig.  18  Stability  of  Eqs  ( 28-29), Or. 1 
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Fig.  19  Stability  of  Eqs. (28-29), Or, 2 
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i slo’je  »=  -l/N 
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x « 0 
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,-4''  7=0 

y « 0 i = 0 

Fig.  20  Isoclines  for  Eqs.  ( 20-29 ) 

M 
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saddle 


saddle 
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-i I 
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stable  focus  or  node  for  R > 0 
Fig.  21  Stability  of  Eqs. (23-29),  Or.  3 


G » R/S 


i 


A 


k - f (y  ) , f (x  ) = k 

x xVJs  ■ yv  s'  y 


A - -(a^)ysfy'(xs)  “ -SM 


B = 0 

C = -(cV'VXs  = ~R 
D - " ”RM 

where  R = ci^x^/  G «*  R/S 

s - v A 

M ■ y <*3>  » - y (y.) 

The  isoclines  for  y = 0 in  this  case  are  y - 0 and  f (x  ) = k , 

J y - t>  y 

and  for  x » 0 are  x = 0 and  x = kx  - f (y) , These  isoclines  are 
plotted  in  Fig,  20,  where  fx(y)  has  been  chosen  to  give  a parabolic 
curve  and  f (x)  has  been  chosen  to  give  a pair  of  vertical  lines. 

Four  singularities  belonging  to  Gr,  3 exist.  Quantity  N can  be  found 
from  the  slope  of  the  isocline  for  x = 0 at  the  singular  point,  as 
was  described  in  discussing  the  most  general  case.  Quantity  H must 
be  calculated  from  its  defining  equation. 

The  stability  diagram  for  Gr.  3 is  shown  in  Fig,  21. 


V.  Degenerate  Cases 

One  degenerate  case  is  that  described  by  the  equations 

x = (a:/kx)(kx  “ x - y)x  (34) 

y « (cy7 ky)(ky  - y - x)y.  (35) 

These  equations  are  of  the  most  general  type,  but  the  coefficients  are 

particularly  simple.  The  isoclines  for  y = 0 and  x ~ 0,  other  than 

the  two  axes,  are  straight  lines  parallel  to  one  another,  so  that  there 

are  no  singularities  in  Gr.  4.  If  > 0,  a > 0,  and  k^  > ky  > 0. 

the  only  stable  singularity  is  that  at  iV.  , “ k , y = 0,  and  the 

S 2C  s 

solution  curves  are  as  shown  in  Fig.  22. 
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A second  degenerate  case  occurs  if  a 


32. 

k “ k 


a = a and  k = __ 
x y x y 

in  Eqs.  (34)  and  (35).  This  is  the  only  combination  of  parameters  in 
the  general  equations  for  which  [kx  - x - f^XyJj  s fk^  - y - f^(x)J  sc 

that  dy/dx  3 y/x.  If  a > 0,  the  origin  is  an  unstable  node.  An 
infinity  of  stable  points  are  located  along  the  line  y = k - x,  and 
the  solution  curves  are  shown  in  Fig.  23. 


VI.  General  Effects  of  Parameters 

Some  comments  can  be  made  about  the  effects  of  the  parameters  that 

appear  in  the  general  equations,  Eqs.  (l)  and  (2).  The  coefficients 

ax  and  a are  the  basic  quantities  in  determining  the  growth  rates. 

Large  positive  values  of  these  quantities  tend  to  produce  rapid  growth. 

Quantities  and  k determine  the  final  values  of  x and  y 

when  a and  a are  positive  and  there  is  no  coupling  between  x 
**  7 

and  y. 

The  functions  f (y)  and  f (x)  determine  the  shape  of  isocline 

x y 

curves  for  x - 0 and  y - 0,  respectively.  If  these  functions  are 
those  of  Eqs.  (3)  and  (4),  the  isoclines  are  parabolic  in  shape.  The 
parabola  for  x *=  0 always  passes  through  the  point  x = k^,  y ~ CL 
At  this  point  the  slope  of  the  parabola  is  dy/dx  = -l/p  . If  p 
varies,  the  parabola  moves  as  shewn  in  Fig.  24,  but  its  shape  does  not 
change.  The  shape  is  determined  by  y?,  and  the  effect  of  changes 
in  yx  is  shown  in  Fig.  25.  The  dotted  line  in  each  of  these  figures 
is  the  locus  of  the  vertex  of  the  parabola  as  it  moves. 

One  type  of  solution  that  is  of  particular  interest  is  that  which 
is  oscillatory  in  nature.  The  simplest  pair  of  equations  that  gives 
oscillations  is  the  following 
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Fig.  ZU  Parabola  Cor  x - 0,  p varies 


Fig.  25  Parabola  fox-  - 0,  y varies 
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X « (a  - by)x 


y " -Cay  - byx)y . (37) 

These  are  examples  of  Sqs,  (18)  and  (19),  where  a * a , a •=  -a  , 

x x y y 

Oj/kjj.  - bx>  3y/ky  - -by#  fx(y)  - y,  fy(x)  = x.  A singularity  exists 

at  x„  ” a /b  = k , y »'  a /b„  » k , For  small  changes  near  this 
s y y y*  " s x x x 

singularity,  the  following  equations  apply 

u - -(a^b^/b^v  (38) 

v = (axby/bx)u  (39) 


to  that 


u + a c u = 0.  (40) 

x y 

The  solution  for  this  equation  is  a periodic  oscillation  having  the 

angular  frequency  u>  = (a  a ) ' . In  terms  of  the  coefficients  used 

x y 

originally,  a and  a determine  the  frequency  of  oscillation  while 

J\.  (V 

k and  k locate  the  mean  values  of  x and  y about  which  oscil- 
x y 

lation  takes  place. 

Only  linear  terms  appear  in  the  functions  and  fy  in  Eqs.  (36) 

and  (37).  If  additional  tenns  are  introduced,  a more  complicated  solu- 
tion results,  A simple  example  illustrating  the  effects  has  extra 
terms  in  only  one  of  the  two  equations, 

X = (ax  - b^y  - c,y  - dxx.)x  (41) 

y - -(ay  - byx)y . (42) 

The  singularity  about  which  oscillation  may  occur  is  located  at  x , y , 

S 5 


given  by 


a - bv  - c y - d x = 0 
x rs  xr  s x s 


a -bx.  =0;  x * a /b  , 
y y s ’ s y y 

For  small  changes  near  this  singularity 

u = -d  x u - (b  x + 2c  x y )v 

X S X S X SJS 

v b y u 

r s 


L. 
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so  that 


~r 


35, 

*ii  - dxxsu  - ay(ax  ♦ V-32  - dxxs)u  - 0.  (47) 

2 

The  term  dxx  in  Ec,  (41)  introduces  damping,  while  both  this  term 

p 

and  the  term  c^y  ~x  change  the  location  of  the  singularity  and  modify 
the  frequency  of  oscillation. 

In  general,  quadratic  terms  in  functions  f and  f , as  in 

x y 

EqSt  (3)  and  (4),  lead  to  a modification  of  the  frequency  of  oscil- 
lation, to  a change  in  the  location  of  the  singularity,  and  also, 

2 * 

perhaps,  to  additional  singularities.  So  long  as  the  terms  x and  y 

are  missing,  as  in  Eqs,  (1C)  and  (19),  coefficients  B arid  C are 

identically  zero  near  the  singularity  whsre  oscillation  occurs.  Thus, 

there  can  be  no  damping  and  the  oscillation  remains  periodic. 

In  the  more  general  case  of  Eqs.  (l)  and  (2),  oscillation  may 

occur  also,  as  was  true  in  Figs.  11  and  12.  However,  the  presence  of 
2 2 

the  terms  x and  y introduces  damping,  the  algebraic  sign  of  which 
depends  upon  the  signs  and  magnitudes  of  terms  in  the  equations.  In 
general,  quantity  (B  + C)  is  not  zero,  and  the  oscillation  either 
builds  up  or  decays,  bith  a particular  adjustment  of  coefficients,  it 
is  possible  to  make  (B  + C)  « 0,  in  which  case  a periodic  solution 
would  exist.  Any  small  changes  in  coefficients  would  destroy  this 
periodicity,  however.  The  amplitude  of  the  peric.  ic  solution  would 
depend  upon  initial  conditions. 

It  does  not  appear  that  a limit  cycle,  representing  a periodic 
solution  with  amplitude  determined  solely  by  the  equations,  can  exist, 

A limit  cycle  might  occur  about  an  unstable  focal  point  if  a suitable 
positive  damping  effect  came  into  play  as  the  amplitude  of  the  oscil- 
lation increased.  In  all  the  examples  that  have  been  studied,  an 


I 

36. 

unstable  focal  point  always  leads  to  a solution  that  ultimately  runs 
away,  as  in  Fig.  11. 

All  of  the  preceding  discussion  has  been  con-,  wind  with  the 
solution  of  Eqs.  (l)  and  (2)  in  which  only  the  relation  between  variables 
x and  y at  given  instants  in  time  is  considered.  Curves  of  y as  a 
function  of  x have  been  obtained  with  time  as  a parameter,  but  a 
scale  of  time  along  these  curves  is  not  available.  In  general,  it 
appears  to  be  quite  difficult  to  find  a solution  for  x,  say,  as  a 
function  of  time.  Such  a solution  would  require  that  y be  eliminated 
from  the  two  equations,  leaving  a single  equation  in  x and  t»  This 
single  equation  would  be  of  second  order  and  would  contain  a number  of 
nonlinear  terms.  In  all  but  very  simple  cases  such  as  Eqs.  (36-3?) 
and  (Zil-42),  the  solution  of  this  equation  by  conventional  analytical 
methods  appears  to  be  a hopeless  task. 
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